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Abstract 

^ Recent data from long-baseline neutrino oscillation experiments have provided new information 
O-i on 6*13, hinting that 0.01 < sin^2^i3 < 0.1 at 2 a confidence level. In the near future, further 
(~| confirmation of this result with high significance will have a crucial impact on the optimization 
of the future long-baseline neutrino oscillation experiments designed to probe the neutrino mass 
^ ordering and leptonic CP violation. In this context, we expound in detail the physics reach of an 
\^ experimental setup where neutrinos produced in a conventional wide-band beam facility at CERN 
^ are observed in a proposed Giant Liquid Argon detector at the Pyhasalmi mine, at a distance 
\^ of 2290 km. Due to the strong matter effects and the high detection efficiency at both the first 
Q\ and second oscillation maxima, this particular setup would have unprecedented sensitivity to the 
O neutrino mass ordering and leptonic CP violation in the light of the emerging hints of large ^13. 
^ With a 10 to 20 kt 'pilot' detector and just a few years of neutrino beam running, the neutrino 
j> mass hierarchy could be determined, irrespective of the true values of 5cp and the mass hierarchy, 
*^ at 3(T (5cr) confidence level if sin^ 26*13 (true) = 0.05 (0.1). With the same exposure, we start 
j-j to have 3o" sensitivity to CP violation if sin^ 2^13 (true) > 0.05, in particular testing maximally 
CP-violating scenarios at a high confidence level. After optimizing the neutrino and anti-neutrino 
running fractions, we study the performance of the setup as a function of the exposure, identifying 
three milestones to have roughly 30%, 50% and 70% coverage in Scpitine) for 3 a CP violation 
discovery. For comparison, we also study the CERN to Slanic baseline of 1540 km. This work 
nicely demonstrates that an incremental program, staged in terms of the exposure, can achieve 
the desired physics goals within a realistically feasible timescale, and produce significant new 
results at each stage. 
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1 Introduction and motivation 



Over the last thirteen years or so, marvellous data from world-class experiments involving neutrinos 
from the sun U\ff\, the Earth's atmosphere [SlIO], nuclear reactors 10 11 , and accelerators 12 13 



have firmly established the phenomenon of neutrino oscillations. This implies that neutrinos have 
mass, an exclusive example of experimental evidence for physics beyond the Standard Model. The 
T2K experiment in Japan is a new addition to this remarkable list of successful experiments and 
has recently provided an important breakthrough in the 6'13-driven z/^ — )► u^. appearance oscillation 
channel. 

The T2K experiment has reported the first indication of electron neutrino appearance in the 
Super-Kamiokande detector, located at a distance of 295 km from an accelerator-based i/^ source. 
The neutrinos have an energy of roughly 0.6 GeV and are produced at J-PARC 14 . The obser- 



vation of six electron-like events with an estimated background of 1.5 events is consistent with 
subdominant z/^ — )■ z/g transitions in the atmospheric sector and indicates a non-zero value of ^13 
at 2.5 cr significance. Within a couple of weeks of the T2K results, the MINOS collaboration an- 



nounced the observation of 62 electron-like events with an estimated background of 49 events 15 
which excludes a zero ^13 at 8 



i C.L. (1.7(t). 

These new results have now been included in global fits of all available neutrino oscillation 
which point towards a non-zero value of 613 at around 3 a C.L. and hint that sin^ 2^13 



16 



17 



data 

lies between 0.01 and 0.1 at 2cr C.L. with a best-fit value of 0.051 (0.063) for normal (inverted) 



mass ordering 17 



With further running of the T2K |18 and MINOS 15 experiments, and the impressive start 



up of reactor neutrino experiments (Double Chooz 19 , RENO 20 , and Daya Bay ^Tj) and of the 



upcoming accelerator experiment NOz/A [22], it is likely that the positive evidence for ^13 will be 
confirmed with high statistical significance in the next few years. A non-zero, and in particular a 
large value of ^13 will provide a boost to the field of neutrino oscillation physics, making it possible 
to discover leptonic CP violation (CPV) if the Dirac CP phase, 6cp, is not equal to 0° or 180°, 
and allowing for a direct determination of sgn{Amlj) aka, the neutrino mass hierarch}|^ 

Given their relatively short baselines, narrow band beams and limited statistics, the present 
generation T2K and NOz/A experiments have a limited reach in probing CPV and the neutrino 



mass hierarchy ^23j, even for the large values of 613 indicated by the recent global fits 16,17 



Hence, new long-baseline experiments with intense neutrino beam sources and advanced detector 



technologies are needed 24 to fathom the hitherto uncharted parameter space of the neutrino 



mixing matrix well beyond the capabilities of T2K and NOz/A. 

the ability of future long-baseline neutrino 



As discussed in the literature (see e.g. Refs. 24 25] 
experiments to discover CPV and the mass hierarchy depends on the achievable event statistics 
and hence strongly on the value of ^13. The definitive confirmation of large ^13 will have a 
significant impact on the realization of future long-baseline neutrino oscillation experiments, for 



^In this paper we define Am?- 



and refer to sgn{Am^i) as tlie neutrino mass hierarclry - sgn^Arri^i) > 



is called 'normal hierarchy' (NH) while sgn{Am^i) < is called 'inverted hierarchy' (IH). The arguments are 
valid for both hierarchical as well as quasi-degenerate neutrino mass spectra. 
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which planning has up to now focused on a staged approach to achieve sensitivity to increasingly 
smaller values of 61^. This approach has been exemplified in the optimization of the Neutrino 



Factory 26 and Beta-beam 27,28 for which it is possible to discover a value of sin^ 26'i3 as small 
as 10~^. However, following the evidence for a relatively large value of 613, the focus of future 
optimizations will be the possibility to explore CPV and the mass hierarchy for a given value of 
sin^ 2^13. 

Accordingly, in this paper, we analyze the impact of large ^13 on the optimization of future long- 



baseline superbeam experiments with a special emphasis on potential European scenarios 29 32 
The technology of conventional beam is well established and understood, and when capable of 
sustaining very high proton beam powers (above a megawatt), they are known as superbeams. 
Superbeam experiments are 'super' in the sense that they use proton beams of unprecedented 
strength and detectors with huge fiducial masses (tens or hundreds of kt). In a superbeam, we 
shoot a target with as many protons as can be provided, and then the produced mesons (mostly 
pions and some kaons) are focused and sign-selected using a magnetic horn. Henceforth, we will 
denote the beam produced by positive horn focusing as 'PHF' and the beam created by negative 
horn focusing as 'NHF'. In our study, we rely on a wide-band beam (WBB) approach, a type 
of superbeam where most of the produced pions are utilized and the on-axis detector is exposed 
to a wide range of energies. This results in higher statistics in comparison with an off-axis or a 
narrow-band beam and also allows several L/E values to be studied in a single experiment, which 
helps to resolve parameter degeneracies |33]. The main caveat with the WBB approach is the 
presence of a high-energy tail in the beam which enhances the feed-down of neutral current (NC) 
background rates. Therefore, an advanced detector technology with good energy resolution and 
good NC rejection capability is vital to use the broad spectrum efficiently. 

In this paper, we focus on the CERN to Pyhasalmi baseline of 2290 km[^ which is the longest 



baseline considered in the LAGUNA design study 31,32 . An outcome of this study is that the 
Pyhasalmi mine in Finland has the appropriate infrastructure, at the right depth, to house an 
underground 100 kt liquid argon detector [36j. The baseline is also sufficiently long such that 
matter effects are quite strong and we can expect that the mass hierarchy can be determined 
at a high confidence level if ^13 is large, irrespective of the true value^ of the other oscillation 
parameters, much earlier than before we start to have sensitivity to CPV. Additionally, the first 
and second oscillation maxima for the CERN-Pyhasalmi baseline are at 4.5 GeV and 1.5 GeV for 
Am|^ = 2.45 X 10~^eV^, well above the typical threshold of a detectoi]^ This enables us to 
observe the first and second oscillation maxima with high efficiency, and also the uncertainties 
on cross-sections due to nuclear effects are negligible at these energies. Higher energies are also 
preferred experimentally as interactions are easier to identify and the total cross-section increases 



34 
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^Note that this baseline is very close to the 'bimagic' distance 

■^We distinguish between the 'true' values of the oscillation parameters, which are the values chosen by Nature 
(or the values we assume that Nature has chosen and which we use to compute the data), and their fitted values. 
Throughout this paper we indicate the true value of a parameter by writing '(true)' after that parameter. 

^Note that the first oscillation maximum is also above the energy of the tau kinematic threshold. This class of 



events was studied in Ref. 37 . 
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with energy. 

We compare the physics reach of the CERN to Pyhasalmi basehne with the CERN to Slanic|^ 
basehne of 1540 km, the second longest basehne option in Europe. This basehne is of a fairly 
similar distance to the Fermilab to Homestake baseline of 1290 km, being actively studied for the 
Long-Baseline Neutrino Experiment (LBNE) [33||38}|39] in the USA. We present a comparison 
between the CERN-Pyhasalmi and CERN-Slanic baselines, in terms of sensitivity to the mass 
hierarchy and CPV. 

As far as the choice of far detectoi[^ is concerned, here we will focus on a Giant Liquid Argon 



Charge Imaging Experiment with a target mass up to 100 kt, as proposed in Refs. 40,41 . An 



initial 'pilot' version could have a mass of 10 to 20 kt. The Liquid Argon Time Projection Chamber 



(LArTPC) technique 42 offers an excellent environment in which to reconstruct electron events in 
the GeV range and above with high efficiency, while considerably suppressing the NC background. 

A relatively large value of allows us to pursue a staged approach in terms of the size of 
the experiment. For a given experimental setup, the physics sensitivity is mainly dictated by the 



incident proton beam power on the target 43 , and by the fiducial mass of the far detector. In this 
study, we assume a CERN-based scenario with a primary proton energy of 50 GeV[^ |30] which 
delivers 3 x 10^^ protons on target (pot) per year, using a beam power of 1.6 MW and running 
for 200 days per calendar year. In order to demonstrate our staged experimental strategy, we 
introduce the quantity 'exposure' which is defined as 

exposure = (potperyear) x (fiducial mass of detector in kt) x (total runtime in years), (1) 

and has units of pot-kt. The maximum detector mass of 100 kt combined with a total ten years 
of running (combined PHF and NHF) gives a total exposure of 3000 x 10^^ pot-kt. We treat the 
exposure as a free parameter in our study with an upper limit of 3 x 10^'* pot-kt and quantify the 
performance of both baselines in terms of the exposure. The following are the main topics that 
we address in our study, in the light of the emerging hints of large ^13: 

• Since the determination of sgn{/\m\^ can be accomplished with a smaller exposure than 
that which is needed for substantial CPV discovery, we first optimize the sharing of the total 
exposure between PHF and NHF in order to gain a 3 a and 5 a confirmation of the mass 
hierarchy, for both choices of the true hierarchy and for 100% of the values of 5c'p(true). 
The information obtained about sgn{^m\i) at the early stage of the experiment plays a 
vital role in maximizing the sensitivity of the setup to CPV - the right choice of the mass 
hierarchy is essentially a mandatory input for discovering CPV. Possibilities to determine it 



with atmospheric neutrinos [44 ,45 have been explored, but are less efficient and more prone 



to systematic errors than the direct method considered here. 



■''The Slanic site in Romania provides in practice only a shallow depth, which would greatly limit the sensitivity 
to astrophysical sources and proton decay. 

^In the real experiment, we will also have a small near detector to reduce the impact of systematic uncertainties. 
''We plan to perform a similar study with a primary proton beam energy of 400 GeV. 
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We propose an incremental approach by identifying milestones of exposure which will enable 
the value of Sqp to be gradually constrained, by incrementally increasing the CPV discovery 
coverage from 30% to 50% and finally 70% of the values of 5cp(true), at a confidence level of 
3 a. We believe that this approach is fully justified by the fact that at present, the favorable 



values (e.g. close to maximal) of 6cp{trne) are still allowed 16,17 . If Nature is as kind as it 
appears to be for 6^13 such that 5c'p(true) is not close to CP conserving values 0° or 180°, then 
we will be able to discover the CPV quite early in the time frame. An incremental approach 
is also justified in view of the challenges (some of them are unknown) involved in operating 
very high power superbeams and in building giant underground neutrino detectors, which 
makes such an approach effectively safer and possibly more cost-effective. 

The paper is organized as follows. We begin with the relevant oscillation probabilities in Section 
2. In Section 3, we give the details of the experimental setup. We briefly describe the CERN-based 
superbeam facility and the fluxes that we use for the Pyhasalmi and Slanic baselines. Then we 
discuss the characteristics of the LArTPC. In Section 4 we present the expected event rates for 
the CERN-Pyhasalmi and CERN-Slanic baselines. The details of our numerical technique and 
analysis procedure are discussed in Section 5. In Section 6, we present our results. Finally, we 
summarize and draw our conclusions in Section 7. 



2 The P^e oscillation channel 

The use of an appearance channel in which the neutrino changes flavor between production and 
detection is necessary in order to probe CPV in neutrino oscillations. In a superbeam, we deal 
with the transition probabilities P{iy^ — ?■ t'e) and -P(t'/i — J- i^e) in order to explore the neutrino mass 
ordering and CPV. For both the baselines considered in this work, matter effects are significant 
and the exact expressions of the three-flavor oscillation probabilities including matter effects are 
very complicated. Therefore, to demonstrate the nature of neutrino oscillations as a function of 
baseline and/or neutrino energy, it is useful to have an approximate analytic expression for P^e 
(the T-conjugate of Pep) in matter 



quantities 6'i3 and a = Aml^/ Aml^ 
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keeping terms only up to second order in the small 
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where 



A^^^, A = ±2V2G,N,E. (3) 



Eq. (|2j) has been derived under the constant matter density approximation. Here, A is the matter 
potential, expressed in terms of the electron density, Ne-, and the ( ant i-) neutrino energy E and is 
positive for neutrinos and negative for anti-neutrinos. The ' — ' sign which precedes the term C_ 
refers to neutrinos whereas the '+' refers to to ant i- neutrinos. The first term of Eq. ([2]) (Cq) can be 
used to obtain information about the value of ^13. This is also the term which contains the largest 
Earth matter effect and which can therefore be used to measure the sign of Amg^^. The term Ci 
is independent of both ^13 and 5c p and depends mainly on the solar parameters, Am2i and 9i2- 
The term C_ is the CP- violating part. The term C+, although ^^p-dependent, is CP-conserving. 
While we will use this formula to discuss our results in some cases, all the simulations presented 
in this paper are based on the full three-flavor neutrino oscillation probabilities in matter, using 
the Preliminary Reference Earth Model (PREM) [Sl] for the Earth matter density. 

In spite of using both neutrinos and anti- neutrinos, a serious complication for all long-baseline 
experiments involving these channels emerges from discrete degeneracies which appear in three 
forms: the (6'i3,5cp) intrinsic degeneracy ^21, the {sgn{^'m\i),5c^) degeneracy [53], and the 



(^23, 7r/2 — 6*23) degeneracy 54 . This leads to an eight-fold degeneracy [55], with several degenerate 
solutions in addition to the true one. The presence of these degenerate solutions can severely 
reduce the sensitivity of an experiment. For both the 2290 km and 1540 km baselines, we can 
extract information from the first as well as the second oscillation maximum quite efficiently using 
the LArTPC detector, and this helps to evade the problem of intrinsic 52 and sgn{^m\i) [53] 



ambiguities which can cause the vr-transit 56 effect even for large values of sin 2^13 (true). 



3 Description of experimental setup 

In this Section, we first present the CERN-based superbeam fiuxes. Then we describe various 
properties of the LArTPC detector which are relevant to this work and summarize them in Table [Tj 



3.1 CERN-based superbeam fluxes 

A conventional (dominantly muon) neutrino beam is produced from the decay of secondary hadrons 
(pions, kaons) which are created by high-energy primary protons impinging on a long target. The 
neutrino beam direction is defined by the orientation of the decay tunnel downstream of the 
target. An important feature is the insertion of focusing elements (so called horns), whose shape 
and current density can be chosen to direct emerging secondary hadrons of a given momentum 
into the decay tunnel, thereby defining the energy spectrum of the resulting muon neutrino beam. 
The magnetic property of the horn is used to select the sign of the produced mesons. When the 
horn acts in the positive (negative) mode, then the resulting beam is known as 'PHF' ('NHF'). 



5 



o 
o 



10° 



Q. 



g_10' 

CO 

cu 

X 
3 



10'' 



v^ 2290 km (PHF) 




5 6 7 8 9 10 



Energy [GeV] 



lO'' 



o 
o 



.10° 



cu 

Q. 



g_10' 

CO 
CU 

>. 

X 
13 



10" 




3 4 5 6 7 
Energy [GeV] 



9 10 



10" 



o 
o 



10° 



CU 



g_10' 

CO 

cu 
>. 

ho' 



lO"" 



2290 km (NHF) 




1 osc max 



1 23456789 10 



Energy [GeV] 




1 23456789 10 
Energy [GeV] 



Figure 1: The upper panels depict the un-oscillated PHF (left panel) and NHF (right panel) superbeam flux 
spectrum arriving at a detector placed at Pyhasalmi (a distance of 2290 km from CERN). The lower panels show 
the same for the CERN-Slanic baseline (a distance of 1540 km). The blue dot-dashed and the orange dotted vertical 
lines display the locations of the first and second oscillation maxima. 
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The customary approach for long-basehne oscillation experiments is to tune the muon neutrino 
flux around the energy of the first oscillation maximum for a given baseline. For our studies, which 
give priority to the measurement of the oscillatory pattern of the flavor conversion processes as a 
function of the neutrino energy, a wide-band neutrino beam tuned at the first oscillation maximum 
and broad enough to provide coverage of the second maximum is the adequate solution. 



Presently the CERN to Gran Sasso Laboratory neutrino beam (CNGS) [57,58 is produced 
using 400 GeV protons from the SPS accelerator and the neutrino flux is peaked around 17 GeV 
to optimize the rate of CC events in the far detector. 

Let us now consider options for the CERN-Pyhasalmi and CERN-Slanic baselines. Prospects 
for high-intensity proton sources at CERN and their relevance for high intensity neutrino beams 



were discussed in Ref. 30 . The existing 400 GeV SPS accelerator, profiting from the LHC lumi- 
nosity upgrades, will be able to deliver beams with higher power in the near future, but reaching 
megawatt beam power at the SPS will require significant replacement of its injector chains and 
a stretch of its performance. A more appropriate solution which was advocated is to consider a 
new high power accelerator (HP-PS2), with a proton energy in the range 30 to 50 GeV |30j, with 
for example an intensity of 2.5 x 10^^ protons and a fast cycle of 1.2 seconds, corresponding to 
a proton beam power of 1.6 MW. A conceptual design is being developed in the context of the 



LAGUNA-LBNO design study 31 , 32 



In this paper, we have assumed that a 50 GeV HP-PS2 drives the neutrino beam for 200 days 
per calendar year, corresponding to an integrated 3 x 10^^ pot per year. The neutrino fluxes were 



computed by A. Longhin 59 , 60 . Also, these fluxes have been optimized for each baseline to have 
a peak at the first oscillation maximum as can be seen from Fig. [1} The blue dot-dashed vertical 
line in each panel of Fig. [T] shows the location of first oscillation maximum which is at 3 GeV 
(4.5 GeV) for the CERN-Slanic (CERN-Pyhasalmi) baseline for Am^^ = 2.45 x 10"^ eV^ The 
left panels of Fig. [l] display the flux spectrum produced by PHF where is the dominant beam 
component and the right panels show the flux spectrum obtained from NHF where dominates. 

The superbeam fluxes suffer from a hurdle of an intrinsic (^e) contamination, which poses 
a serious problem in the form of backgrounds in — v^. {i?^ — t- z/g) appearance study. At the flux 
level, we have 0.62% (0.45%) intrinsic contamination of z/g (^'e) compared to z/^ (z/^) in the beam 
for CERN-Pyhasalmi baseline. In the case of CERN-Slanic baseline, the z/g (z/g) contamination is 
at the level of 0.57% (0.39%) compared to (z/^) in the beam. 

3.2 The LArTPC detector 

Over the last several years, dedicated R&D has been pursued to have a new giant next-generation 
and multi-purpose neutrino observatory based on the LArTPC technology with a total mass in 



the range of 100 kt 40 ,41 ,61 -65 , devoted to particle and astroparticle physics, providing new and 
unique scientific opportunities in this field and likely leading to fundamental discoveries. 

The LArTPC is a powerful neutrino detector for uniform and high accuracy imaging of massive 



active volumes 42 . It is based on the fact that in highly pure Argon, ionization tracks can be 



drifted over distances of the order of meters. Imaging is provided by position-segmented electrodes 
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Detector characteristics 


LArTPC (Both /i± and e^) 


Fiducial mass 


10 - 100 kt 


Neutrino energy threshold 


0.5 GeV 


Detection efficiency (e) 


100% for /i± 
80% for e± 


Energy resolution (SE) (GeV) 


0.15v/E/GeV for CC /i± and e± sample 


NC background smearing 


Migration matrices (different for u and u) 


Bin size 


0.125 GeV 


NC background rejection efficiency 


99.5% 


Background from misidentified muons 


0.5% 


Efficiency for intrinsic t'e/z^e contamination 


80% 


Signal error (systematic) 


5% 


Background error (systematic) 


5% 



Table 1: Detector characteristics used in the simulations. The last two rows show the reference 



choice of systematic errors on signal and background. Details can be found in Section 3.2 



at the end of the drift path, continuously recording the signals induced. Both single and double- 



phase LArTPCs 66 have been successfully operated. Reference event timing (Tq) is provided by 
the prompt scintillation light. Application of this technology, originally developed at CERN, to 
large detectors was pioneered by the ICARUS effort which culminated in the successful operation 



of the T600 on surface 42 . Underground operation was recently reported in Ref. 67 . Motivated 
by the necessity of massive detectors for neutrino physics and proton decay search, several designs 
of large LAr detectors (much more massive than ICARUS T600) have emerged in the literature 



during the last decade 40 , 41 , 68 - 72 



In our study, the target mass of the detector is not fixed and it starts from a 'pilot' of 10 kt 
and can go up to 100 kt. We consider the reconstructed neutrino energy threshold to be 500 MeV|^ 
for all the channels. We use 80% detection efficiency, e, for charged current (CC) electron or 
positron events (see Table [T]). For /i^ events, the efficiency is 100%. The energy resolution 
of the detector is expected to be very good and we assume that the neutrino energy will be 
reconstructed with a GauBian energy resolution of 15% of ^^/E/GeV. We have considered the Ue 
and the z/g appearance channels in our study, with the background contribution coming from the 
intrinsic i^e/'^e contamination of the beam, the number of muon events which will be misidentified 
as electron events, and NC events. We assume a rejection efficiency of 99.5% for NC eventajand 



®In a LArTPC, one can actually detect energies < 100 keV, but at these energies the triggering becomes 
problematic, even in an underground location. Also, at lower energies the cross-sections are more uncertain due to 
nuclear effects. 

^Our choice is conservative in the sense that the excellent spatial resolution and the capability of observing very 
low-energy particles might enable a LArTPC to reject NC events even more efficiently. 
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the correlation between the visible energy of background NC events and the neutrino energy is 



implemented by migration matrices 73 . To estimate the backgrounds coming from the intrinsic 
z/g/t'e contamination, we use the same 80% efficiency that we have used for our signal. We also 
assume that there is a 0.5% chance of misidentifying muons as electron events. We also add the 
information coming from the v^/v^ disappearance channels. For these channels, NC events are 



the main source of backgrounc ^ We do not consider backgrounds due to atmospheric neutrinos 
as the energy of the beam events is well above that of the an average atmospheric event and also, 
the timing information and directionality of the beam are sufficient to reduce these backgrounds 
to an insignificant level. We also include signal and background normalization errors (both 5%), 
uncorrelated among all oscillation channels. 



4 Event rates for the 2290 km and 1540 km baselines 

In this section, we present the expected event rates in the LArTPC for the CERN-Pyhasalmi 
and CERN-Slanic baselines using a reference exposure of 1500 X 10^1 pot-kt (see Eq. ^) for each 
channel. 

The number of electron (positron) event^^ in the i-th energy bin in the detector is given by 
^i = ^rij / dEAHE)a,SE)R{E,EA)P,e{E), (4) 

4:71 L Jq 

where T is the total running time, is the number of target nucleons in the detector, e is the 
detector efficiency and R{E, Ea) is the GauBian energy resolution function of the detector. For 
electron (positron) events, a^^ is the neutrino (anti-neutrino) interaction cross-section. These 



cross-sections have been taken from Refs. 74,75] and they are estimated for water and isoscalar 
targets. However, we will need to correctly evaluate the cross-sections for argon when we analyze 
the real experiment. The quantities E and Ea are the true and reconstructed (anti-)neutrino 
energies respectively, and L is the baseline. 

Fig. [2] portrays the expected signal and background event spectra as a function of reconstructed 
neutrino energy including the efficiency and background rejection capabilities for the 2290 km 
(upper panels) and 1540km (lower panels) baselines, using an exposure of 1500x10^^ pot-kt in 
each channel. Results are presented for sin^2^i3 = 0.05 and (5cp = 0° with a NH; for all other 
oscillation parameters, we have used the benchmark true values as given in Table |4j The left 
panels show the performance of the Ve appearance channel and the right ones are for the Ve 
appearance channel. In all the panels of Fig. |2| one can clearly see a systematic downward bias 

"'^'^The [Vf^) 'wrong-sign' contamination of the beam will also act as a background for the (P^) disappearance 
channel. This background is negligible compared to the number of 'right-sign' signal events and moreover, the 
and disappearance oscillation probabilities are similar (except for differences due to the sub-leading matter effect 
terms) which is why they do not cause any problem in this study. 

-^^We also consider events. The number of events can be calculated using Eq. Q, by making appropriate 
changes to the oscillation probability and cross-sections. 
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Figure 2: Expected signal and bacltground event rates including the efficiency and background rejection 
capabilities in the i^e (left panel)/i>e (right panel) appearance channel for sin^ 20i3 = 0.05 and Sep = 0°, as a 
function of the reconstructed neutrino energy. The results are shown for an exposure of 1500 x 10^^ pot-kt for each 
channel. The upper panels are for the CERN-Pyhasalmi baseline and the lower ones for the CERN-Slanic baseline. 
A normal hierarchy has been assumed. For all other oscillation parameters we use the values given in Table |4] In 
all the panels, the blue dot-dashed and the orange dotted vertical lines display the locations of the first and second 
oscillation maxima. The red solid histograms show the total signal+background rates. 
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Channel 


CERN-Pyhasalmi (2290 km) 


CERN-Slanic (1540 km) 


Signal 


Background 


Signal 


Background 


CC 


Int+Mis-id+NC = Total 


CC 


Int+Mis-id+NC = Total 


(NH) 


2364 


419+100+103=622 


3027 


694+210+158=1062 


v^, (IH) 


485 


439+100+103=642 


1110 


710+212+158=1080 


v^, -> j>e (NH) 


304 


128+42+45=215 


622 


192+87+66=345 


Vf, Ve (IH) 


1049 


122+43+45= 210 


1127 


188+89+66=343 



Table 2: Comparison of the total signal and background event rates in the Vf./Vf, appearance channel for CERN- 
Pyhasalmi and CERN-Slanic baselines taking sin^ 20i3 = 0.05 and (5cp = 0°. The results are shown for both the 
choices of the neutrino mass ordering with an exposure of 1500 x 10^^ pot-kt in each channel. Here signal efficiency 
and background rejection factors have been taken into account from Table [ij Here Tnt' means intrinsic beam 
contamination, 'Mis-id' means misidentified muon events and 'NC stands for neutral current. 



Channel 


CERN-Pyhasalmi (2290 km) 


CERN-Slanic (1540 km) 


Signal 


Background 


Signal 


Background 


CC 


NC 


CC 


NC 


v^, (NH) 


19938 


103 


41962 


158 


Vp. Vf, (IH) 


20096 


103 


42486 


158 


i>M ^ '^1^ (NH) 


8488 


45 


17488 


66 


Vf, (IH) 


8545 


45 


17729 


66 



Table 3: Comparison of the total signal and background event rates in the v^/v^j, disappearance channel for 
the CERN-Pyhasalmi and CERN-Slanic baselines taking sin^ 26*13 = 0.05 and (5cp = 0°. For all other oscillation 
parameters we use the values given in Table|4j The results are shown for both choices of the neutrino mass hierarchy 
with an exposure of 1500 x 10^^ pot-kt in each channel. Signal efficiency and background rejection factors have 
been taken into account from Table [ij 

in the reconstructed energy for NC background events due to the final state neutrino included 
via the migration matrices. The blue dot-dashed and the orange dotted vertical lines display the 
locations of the first and second oscillation maxima. The green double-dotted-dashed histogram 
shows the signal event rate and the red solid histogram presents the sum of signal and background 
rate. Although we have enough statistics around the second oscillation maximum for both the 
baselines, its impact is limited due to the fact that the event samples are highly contaminated 
with NC and other backgrounds at lower energies. 

In Table [2| we present a comparison between the total signal and background event rates in the 
z/e/z^e appearance channel for the CERN-Pyhasalmi and CERN-Slanic baselines, taking sin^ 26'i3 = 
0.05 and 5cp = 0°. The results are shown for both choices of the mass hierarchy with an exposure of 
1500x10^^ pot-kt in each channel, including the signal efficiency and background rejection factors 
from Table [H The Earth matter effects in the neutrino channel with a NH drives the number of 
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Figure 3: Event rates in the appearance channel as a function of Sqp- We consider sin^ 26i3 = 0.05. Left panel 
(right panel) portrays the event rates for the CERN-Pyhasalmi (CERN-Slanic) baseline. Here, v (NH) refers to 
the — >■ channel with a NH. Similarly, i? (NH) corresponds to the — > channel with a NH. The results are 
shown for an exposure of 1500x10^^ pot-kt in each case. 



expected events to be very large, compared to the number expected for an IH. Similarly, in the 
anti-neutrino channel, matter effects in the case of an IH mean that the number of predicted events 
is larger than for a NH. At the 2290 km baseline, the impact of matter effects is more significant 
compared to the 1540 km baseline, and therefore the relative difference between the number of 
events for NH and IH is larger for the CERN-Pyhasalmi than the CERN-Slanic baseline. It is clear 
from Table [2] that the most dominant contribution to the background comes from the intrinsic 
Ue/ue beam contamination. Table [3] depicts the total signal and background event rates in the 
z/^/ z/^ disappearance channel for the CERN-Pyhasalmi and CERN-Slanic baselines considering an 
exposure of 1500x10^^ pot-kt in each channel. The information from the v^/v^j^ disappearance 
channel is crucial to further constrain the atmospheric parameters. For this channel, the 1540 km 
baseline has almost twice the number of events that we have in case of 2290 km. 

In Table [2| we show the signal event rates in the Ve/ appearance channel for a specific value 
of (5cp which we have chosen to be one of the CP conserving cases. It is interesting to see how 
the signal event rate varies with the choice of (5cp, which is what we show in Fig. [3j This shows 
the 5cp dependence of the event rates for the CERN-Pyhasalmi (left panel) and CERN-Slanic 
(right panel) baselines. In each panel, the results are shown for both the Vf. and Vf. appearance 
channels considering both choices of the mass hierarchy. Here, we consider sin^ 26'i3 = 0.05 and an 
exposure of 1500x10^^ pot-kt in each case. The number of events varies in a wide range with the 
choice of 5cp for both the baselines. We can see that the maximum (minimum) number of events 
for neutrinos is obtained around -90° (90°) for both NH and IH. For anti-neutrinos, the same is 
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Benchmark Values 


Icr estimated error 


Am2^(true) = 2.45 x lO'^ eV^ (NH) 


a(Am2j = 5% 


Am|i(true) = -2.34 x lO^^ eV^ (IH) 


a(Amii) = 5% 


^23 (true) = 45° 


a(023) = 10% 


Am2^(true) = 7.59 x lO^^ eV^ 


(T(Am2i) = 3% 


012 (true) = 33.96° 


a(^^i2) = 3% 


p(true) = 1 


a{p) = 2% 



Table 4: Chosen benchmark vahies of oscillation parameters and their 1 a estimated errors. In the last row, p 
is the Earth matter density, relative to the value given by the Preliminary Reference Earth Model (PREM). 



true with 5cp — ^cp- We can analyze this behavior with the help of Eq. (|2]). We do not have 
resonant matter effects for either of the baselines and hence A < 1 in Eq. (|3]). It follows that the 
(5cp dependence of the probabihties can be written as P^e — Co + Ci — C_ sin 5cp + C+ cos 5cp 
and Pjie — Cq + Ci + C_sin5cp + C+cos^cp, where Cq.Ci and, C± are independent of 5cp, 
whence the symmetry is manifest. For most of the energies, C_ is positive assuming a NH. Now 
for (5cp = —90° (90°) the term C_sin5cp gives a positive (negative) contribution towards the 
probability for neutrinos. The opposite is true for anti-neutrinos. A similar argument is valid for 
the case of an IH. 



5 Simulation method 



For all the numerical results presented in this study, we have used the GLoBES software 76 , 77 . 



Unless stated otherwise, we have generated our simulated data for the benchmark values given in 
the first column of Table |4} These values correspond to the current best-fit values of the oscillation 
parameters, in light of the current neutrino data 78 . In all fits, these parameters are marginalized 



over, with the 1 a intervals given in the second column of Table |4] included as prior^| (see Eq. (|8 
We also include a 2% uncertainty on the matter density (p) [79[]8l] . 

For our statistical analysis we use the techniques outlined in Refs 
following function: 



56 , 82 and employ the 



Xtotal 



2 _|_ 2 _|_ 2 



The for the z/„ — )■ z/g appearance channel is given by 



prior 



(5) 



i=l 



X, - Xi In + il 



(6) 



^^The setup we are considering can improve upon the current measurements of the atmospheric parameters. 
However, it is useful to add the current informatio to speed up the marginalization. 
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where n is the total number of bins and 



m^}, {^s, ej) = N'^\W}) [1 + vr^e] + N^iiu^}) [1 + 7,%] . (7) 

In the above expression, N^^{{uj}) is the predicted number of CC signal events (calculated using 
Eq. Q) in the i-th. energy bin for a set of oscillation parameters to and N^{{uj}) are the total 
number of background events in bin i. The quantities vr* and vr^ in Eq. ([T]) are the systematic 
errors on the signals and backgrounds respectively. We consider tt* = 5% and vr^ = 5% as 
mentioned in Table [T] The quantities C,s and C,b are the pulls due to the systematic error on the 
signal and background respectively. The data in Eq. ^ enters through the variable Xi = N[^ + N^, 
where iVf^' is the number of observed CC signal events in the detector and A^f is the background, 
as mentioned earlier. We simulate the CC signal event spectrum using Eq. Q for our true values 
of the oscillation parameters given in the first column of Table |4j We have chosen to study a range 



of values of 6^13 such that 0.01 < sin 2^13 < 0.1, consistent with the 2 a range as given in Ref. 17 
We consider the full range of values of 6cp and study both the mass hierarchies. In a similar 
fashion, we calculate the contributions towards Xtotai coming from the other oscillation channels 
(for both neutrino and anti-neutrino modes). In our fit we marginalize over all oscillation 
parameters, as well as the mass hierarchy, by allowing these parameters to vary in the fit and 
picking the smallest value of the function. However, we impose a prior, or external constraint, 
on the atmospheric and solar parameters and on the matter density, through a Xprior^ given by 



x' 



\Aml-^\ — I Am3]^(true)| \ ^ / 623 — 6*23 (true ''^ ^ 



Am|i — Am2i(true)\^ / 612 — ^12 (true) 



+ 
+ 



a(Amii) ; V ^(^12) 



where p is the matter density relative to the PREM density, and the la errors are those given 
in the second column of Table |4j No priors are included for ^13 and 5cp which means that they 
are completely free in the marginalization. However, it is expected that the currently running 
or upcoming accelerator and reactor experiments will not only strengthen the emerging hints of 
non-zero ^13 but also provide an allowed range of sin^ 2^13 at high significance, independent of the 
value of 5cp- This information will be very useful when incorporated in the form of a prior on ^13, 
and will certainly significantly improve the results presented in this study. 

6 Results 

In this section we explore the capability of the facility described in Section |3| with both the 
2290 km and 1540 km baselines, to make a measurement of the true neutrino mass hierarchy (MH) 
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We consider both CC and NC background events; NC events do not depend on the oscillation parameters. 
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and the CPV phase, 5cp (true). For long-basehne experiments, the measurement of the mass 
hierarchy is easier than a measurement of 6cp{tTue) because matter effects enhance the separation 
between the oscillation spectra, and therefore the event rates, of a NH and an IH. Additionally, 
this measurement is one that is 'discrete' as we only need to differentiate between two possibilities. 
A 'discovery' of the mass hierarchy is defined as the ability to exclude any degenerate solution 
for the wrong (fit) hierarchy at a given confidence level; we consider the 3 a and 5 a levels. A 
'discovery' of CPV, if it exists, means being able to exclude the CP-conserving values of 0°, 180° 
at a given confidence level (again, we consider the 3a and 5a confidence levels). Clearly, this 
measurement is increasingly difficult the closer that 6cp{tTue) is to 0°, 180°. Therefore, whilst 
it is possible to discover the mass hierarchy for all possible values of 5cp(true), the same is not 
true for CPV. As such, we focus firstly on the measurement of the mass hierarchy, for which we 
will show that a result can be guaranteed for suitably large ^13, irrespective of the true value of 
(5cp(true), with only a modest exposure. We will then discuss the CPV discovery potential of the 
facility, showing that whilst a larger exposure is needed in order to obtain a reasonable sensitivity 
than for a measurement of the mass hierarchy, it is still feasible to cover a significant portion of 
the parameter space with a moderate exposure. 

We remind the reader that when we refer to 'PHF' (positive horn focusing) and 'NHF' (negative 
horn focusing) in this section, the former refers to a beam consisting predominantly of i/^'s whereas 
the latter refers to a beam consisting mainly of z^^'s, as shown in Fig. [T| 

6.1 Measurement of the mass hierarchy 

We begin by referring back to Fig. [3] where we showed the event rates in the appearance channel 
as a function of 5cp (true), for both a NH and an IH. As an initial simplistic analysis, we can 
firstly look at the relative difference between the number of events for a NH and an IH, for a 
given value of 5cp(true). This gives a rough idea as to how easy it is to distinguish between the 
mass hierarchies, for each of the baselines. We can see that the difference is always greater for the 
2290 km baseline than for the 1540 km baseline, for all values of 6cp{tTue) and both u and z/, which 
is what is expected due to the stronger matter effects at the longer baseline. We also see that the 
distinction between the mass hierarchies is greatest for z/'s around 6cp{tTue) ~ —90° whereas the 
minimum occurs around (5cp(true) ~ +90°. The opposite is true for z/'s, which means that z/'s and 
z/'s have sensitivity to complementary parts of the 5cp(true) parameter space. Therefore, in order 
to minimize the exposure needed to obtain coverage for all values of 5cp (true), it is advantageous 
to combine both PHF and NHF. 

The question is exactly what combination of PHF:NHF is optimal, in the sense of both the 
physics potential and experimental feasibility. Experimentally, PHF is preferable to NHF for 
several reasons. Predominantly, the z/ charged- current cross-section is lower than for u whereas 
the associated background is higher, and the systematic errors for z/ are currently greater than 
for z/. Additionally, the production yield for vr", which decays into z/^, is slightly lower than for 
vr+'s. Therefore we want to know if, firstly, data from NHF is actually necessary, and secondly, 
if it is, the minimum amount of NHF required for optimal performance. In Fig. |4] we show the 
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Figure 4: Total exposure needed to achieve mass hierarchy discovery at 3ct C.L. for 100% values of 6cp{tme) as 
a function of the percentage of NHF in the total exposure, for sin^ 26'i3(true) = 0.05. The results are qualitatively 
similar for other values of sin^ 26*13 (true) and confidence levels. 



total exposure (PHF and NHF combined) needed in order to make a 3 a determination of the mass 
hierarchy for all values of 5c'p(true), as a function of the percentage of the total exposure for which 
the experiment runs in NHF. The results are shown for both the baselines we are considering, and 
for both NH (true) and IH (true). The point at which the lines for NH (true) and IH (true) cross, 
for each baseline, shows the amount of NHF running needed such that the same total exposure is 
adequate for both of the mass hierarchies. Assuming that we have no prior hints as to the nature 
of the true hierarchy before the experiment starts running, then at all other points, the total 
exposure which is needed must be taken to be the higher one. Therefore, without prior knowledge 
about the mass hierarchy, the crossing point is the optimum. 

We find that it is possible to identify the mass hierarchy for all values of Scpitrue) with 
only PHF or only NHF, which indicates that there are no intrinsic degeneracies which cannot be 
resolved simply by accumulating more data. Approximately an order of magnitude more running 
is needed for NHF alone than for PHF alone, in the case of a NH (true), and vice versa for an IH 
(true); a NH favors more data from i/'s whereas an IH prefers data from i/'s. This is a consequence 
of the fact that u oscillations are enhanced by matter in the case of a NH whereas u oscillations 
are enhanced in the case of an IH. 

For 2290 km, the optimum lies at a PHF:NHF ratio of around 30 : 70, whereas the shorter 
1540 km baseline prefers slightly more NHF and requires a higher total exposure. This can be 
understood if we refer back to Fig. [3] which shows that the relative difference between the event 
rates for a NH and an IH is always larger for 2290 km than 1540 km, and hence a lower exposure 
is required for the longer baseline. Additionally, as mentioned earlier, we see how i/'s and u^s have 
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Figure 5: CP fractions for which a mass hierarchy discovery at 3cr (upper panels) and 5cr (lower panels) 
confidence levels is possible as a function of total exposure, with equal sharing of PHF and NHF. Results are 
depicted for three different choices of sin^ 20i3(true): 0.01 (left panel), 0.05 (middle panel) and 0.1 (right panel). 

sensitivity to complementary parts of the Scp{tTue) space. Since we require a minimum threshold 
number of events to perform a 3 a determination of the mass hierarchy, for each value of 5cp (true), 
this means that the 1540 km baseline requires more u data than 2290 km in order to reach the 
threshold number for the values of 6cp{tTue) in which the u data is weakest. 

Bearing in mind that 50% NHF is the maximum that is experimentally favorable, and that the 
difference in exposure required between the exact optimum and the 50 : 50 mixture is small, we set 
a 50 : 50 ratio of PHF and NHF (shown by the dotted vertical black line) to be our reference con- 
figuration for mass hierarchy discovery. In Fig.|4]we show only the results for sin^ 2^i3(true) = 0.05 
at 3 cr confidence but we have verified that the quantitative optimization is practically independent 
of the value of ^13 (true) and the statistical confidence level. 

In Fig. |5] we show how the sensitivity to the mass hierarchy varies as a function of total 
exposure, when we use our reference configuration with equal PHF and NHF. We show the results 
for both baselines and mass hierarchies, and three different values of sin^ 2^13 (true) - 0.01 (left 
panels), 0.05 (middle panels) and 0.1 (right panels). The upper row shows the results at 3cr C.L. 
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Percentage of NHF Percentage of NHF Percentage of NHF 

Figure 6: CPV discovery potential at 3 cr C.L. as a function of the percentage of NHF in the total exposure, for 
sin^ 26'i3(true) — 0.05. Results are shown for three different benchmark choices of total exposure: 300 x 10^^ pot-kt 
(left panel), 1000 x 10^^ pot-kt (middle panel) and 3000 x 10^^ pot-kt (right panel). The results are qualitatively 
similar for other values of sin^ 2013 (true) a-nd confidence levels. 



and the bottom row the resuhs for 5 a C.L.. The dotted vertical black lines mark exposures of 
300 X 10^^ pot-kt and 1000 x 10^^ pot-kt, chosen for reasons that will be explained in the following 
subsection. The right-most edge of the graph lies at our maximum exposure of 3000 x 10^^ 
pot-kt. For all values of sin^ 26'i3(true) there is slightly better sensitivity to a NH than an IH (see 
Fig. |4] where the NH (true) line is at a lower exposure than the IH (true) line for 50% NHF). 
For the longer baseline, sensitivity begins at a lower exposure and the increase in sensitivity as 
a function of exposure is more rapid than for the 1540 km baseline, although the differences are 
not large. Even if sin^ 2^i3(true) ~ 0.01, both baselines can obtain 3 a sensitivity for the complete 
parameter space. The 1540 km baseline is capable of 100% coverage at 5a confidence for values of 
6*13 down to sin^ 26^13 (true) ~ 0.025 (not shown); for the 2290km baseline it is nearly possible at 
sin^ 2^13 (true) ~ 0.01. 

The capability of this setup then, for the large values of ^13 under consideration, is remarkable: 
if we consider sin^ 2^13 (true) = 0.1 and the 2290 km baseline then complete sensitivity at 3cr 
can be gained with an exposure of only 20 x 10^^ pot-kt. This corresponds to, for example, a 
10 kt detector with a beam running for roughly one year. We obtain complete 5 a sensitivity 
with roughly twice this exposure. For sin^ 2^13 (true) = 0.025, the 2290 km baseline is capable of 
achieving 100% sensitivity at 3 cr confidence after roughly 200 x 10^^ pot-kt (20 kt detector running 
for 5 years, for example). 

6.2 CP violation discovery potential 

In Fig. |6]we repeat our optimization of the PHF:NHF running, this time with respect to the CP 
violation discovery potential. Since it is impossible to use 100% coverage as our goal, we instead 
fix the exposure at values of 300 x 10^^ pot-kt (left panel), 1000 x 10^^ pot-kt (middle panel) and 
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3000 X 10^^ pot-kt (right panel) - explained shortly - and plot the CP coverage which is obtained 
when using different ratios of PHF:NHF. Again, we show the results only for sin^ 2^i3(true) = 0.05 
and 3cr C.L. but have verified that they remain qualitatively similar for other values of ^i3(true) 
and confidence levels: altering ^13 (true) effectively just shifts the lines vertically - in the case of 
smaller ^13 (true), the shift is downwards, towards lower CP fractions, and in the case of larger 
^13 (true) the shift is upwards. We see that a true NH prefers more NHF whereas a true IH 
prefers more PHF running, which is the opposite to what we found for mass hierarchy discovery. 
This indicates that whereas the mass hierarchy could reasonably easily be determined simply by 
accumulating sufficient statistics, in the case of CP violation, the most effective measurement is 
made when there are roughly equal numbers of v and v events. This is consistent with the intuitive 
assumption that CP violation is most easily observed by directly comparing z/'s with i/'s. In the 
case of a true NH where the v oscillations are enhanced and v oscillations are suppressed, more 
NHF is required in order to equalize the number of v and v events. The opposite is true in the 
case of a true IH. 

In terms of experimental strategy, this means that one should run with an equal mixture 
of PHF and NHF (shown by the dotted vertical black lines) until the mass hierarchy has been 
determined. Then, in order to achieve the maximal sensitivity to CP violation in the minimum 
possible time, if the hierarchy has been determined to be normal, one should continue running 
with equal PHF and NHF (or with even more NHF if possible) since a NH favors more v data for 
CP sensitivity. However, if the hierarchy has been determined to be inverted, one should decrease 
the amount of NHF to around 30%. This apphes up until an exposure of roughly 1000 x 10^^ 
pot-kt has been reached, after which all admixtures of PHF and NHF with between roughly 20% 
and 80% NHF give similar results, regardless of the mass hierarchy. Therefore, it is possible to 
run with just ~ 20% NHF in this regime and still achieve optimal results. 

Fig. [7] shows the CP sensitivity of the setup as a function of exposure, for sin^ 26'i3(true) = 0.01 
(left panel), 0.05 (middle panel), and 0.1 (right panel), at 3 a C.L. (top row) and 5 a C.L. (bottom 
row). From the central panel, which corresponds roughly to the current best-fit value of sin^ 2^13, 
we define three target exposures from which we can achieve three milestone experimental goals 
of obtaining 3cr sensitivity to CP violation for 30%, 50% and 70% of the (5(7p(true)-space. The 
corresponding exposures are 300 x 10^^ pot-kt, 1000 x 10^^ pot-kt and 3000 x 10^^ pot-kt and are 
marked by the dotted vertical black lines in all the panels. 

In this figure we have considered an equal amount of PHF and NHF, as this is the optimal 
configuration prior to the discovery of the mass hierarchy as explained previously. In general, 
for exposures smaller than the maximum of 3000 x 10^^ pot-kt, the shorter baseline performs 
better, although the longer baseline only requires slightly more exposure (roughly 40% more) 
to attain a particular CP fraction. The difference decreases as the exposure is increased and 
systematic errors become the limiting factor. At the maximum exposure, both baselines perform 
equally well, yielding a maximum CP fraction at 3 a C.L. of between 60% and 75% for both mass 
hierarchies, depending on the value of ^13 (true). At 5(7 C.L. the maximum coverage is 25% to 
60%, again dependent on the exact value of ^13 (true). 

If sin^ 26'i3(true) is just slightly larger than 0.01, the facility begins to have some 3 a sensitivity 
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Figure 7: CP fractions for which a CPV discovery at 3 cr (upper panels) and 5 cr (lower panels) confidence levels 
is possible as a function of total exposure, with equal sharing of PHF and NHF. Results are depicted for three 
different choices of sin^ 26*13 (true): 0-01 (left panel), 0.05 (middle panel) and 0.1 (right panel). 

to CPV even with our lowest targest exposure. If not, and sin^ 2^13 (true) ~ 0.1, our second target 
exposure of 1000 x 10^^ pot-kt is sufficient to produce an appreciable 3 a coverage of ~ 30% for 
both baselines and hierarchies. 

In Fig. [sjwe take our three experimental goals - to cover 30%, 50% and 70% of the 5cp(true)- 
space for CP violation - and ask what the minimum exposure is that is required to achieve these 
goals, as a function of sin^ 26^13 (true). We have again used equal PHF and NHF but the results 
do not change significantly if the ratio is altered slightly. For 30% coverage (left panel), we show 
that our lowest target exposure of 300 x 10^^ pot-kt (the lower dotted horizontal black line), is 
that which can provide 30% coverage in the case that sin^ 26'i3(true) > 0.05 as explained earlier, 
for both baselines and hierarchies. Our next target exposure, 1000 x 10^^ pot-kt (the upper dotted 
horizontal black line), is that which provides 30% coverage if sin^ 2^13 (true) = 0.01. In order to 
reach 50% coverage (middle panel) for this range of sin^ 2^13 (true), we have to consider exposures 
up to 2000 X 10^^ pot-kt (for sin^ 2^13 (true) = 0.01). Reaching 70% (right panel) is feasible with 
our maximum exposure only for certain values of ^^13 around 0.05 - if ^13 is either very small or 
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Figure 8: Total exposure required as a function of sin^ 26'i3(true) to have CPV discovery at 3 a C.L. for 30% 
(left panel), 50% (middle panel), and 70% (right panel) values of (5cp (true). We assume equal sharing of PHF and 
NHF in the total exposure. 



very large, it will not be possible, even with our maximum exposure. Small values of ^13 (true) 
mean that gaining sufficient statistics is a problem, whereas very large values of 6'i3(true) mean 
that the atmospheric term of the oscillation probability is so large that it swamps the CP term, 
effectively acting as a background. However, bearing in mind that the CP coverage which can 



be obtained by a neutrino factory is ~ 80% 26 , the expected performance of this setup is quite 
compelling. 



6.3 Precision measurements of ^13 and Sqp 

Finally, in Fig. [o] we show the fits which can be obtained in the sin^ 2^^i3 — 5cp plane at 68% 
and 95% C.L. (2 d.o.f). We use our three benchmark exposures of 300 x 10^^ pot-kt (left panel), 
1000 X 10^^ pot-kt (middle panel) and 3000 x 10^^ pot-kt (right panel), again with equal sharing of 
PHF and NHF modes. In each plot, we have simulated true values of sin^ 2^13 (true) = 0.01, 0.05 
and 0.1 and 6cp(tYVLe) = 0°, ±90° and 180° (these true values are marked by the black dots), for 
both baselines. Here we show the results for the case of a true NH only; the results for a true IH 
are qualitatively similar. 

The dotted contours which are present for sin^ 2^13 (true) = 0.01 and the two lowest exposures 
correspond to the degenerate solutions which can be obtained by fitting to the wrong hierarchy 
(IH in this case). In other words, for these points, the hierarchy cannot be determined at the 95% 
C.L.. It is important to note that these degenerate regions do not occur at the correct values of 6cp 
(i.e. they do not occur at 5cp(true)), as do the regions corresponding to the correct hierarchy. For 
example, the degenerate solution corresponding to (sin^ 26^13 (true) = 0.01, ^cp (true) = 0°, NH) 
occurs at (sin^ 2^13 (true) = 0.01, 5c'p(true) = —90°, IH) so that the true, CP-conserving solution 
has a degenerate solution at a maximally CP- violating point. Similarly, the degenerate solution 
for the point (5cp(true) = 90° occurs at a CP-conserving point. 
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Exposure = 3000 x 1 pot.kt 

2290 km 68%, 95% C.L. (2 d.o.f.) 

1540 km 

I © ^ 

% ^ ,-rP=^ 

0.05 0.1 
sin^29,3 

Figure 9: 68% and 95% (2 d.o.f.) C.L. contours in the sin^ 26'i3-Jcp plane for sin^ 26li3(true) = 0.01, 0.05 and 
0.1 and (Jcpltruc) = 0°, ±90° and 180° (shown by the black dots). We have simulated a true NH; the results for a 
true IH are qualitatively similar. Equal running of PHF and NHF modes has been assumed with a total exposure 
of 300 X 10^1 pot-kt (left panel), 1000 x lO^^ pot-kt (middle panel) and 3000 x lO^i pot-kt (right panel). 

We find that when statistics are the hmiting factor (small sin^ 2^13 (true) or low exposure), 
the 1540 km baseline produces slightly better results (tighter contours) than the 2290 km baseline, 
because of the higher flux at the shorter baseline. However, once sufficient statistics have been 
accumulated, the longer baseline performs as well as, or sometimes slightly better than, the shorter 
baseline. This is a consequence of the fact that the oscillation peak of the longer baseline covers 
a wider range of energy bins than the oscillation peak of the shorter baseline (recall that we are 
using the same energy bins for the analyses of both the baselines) which in turn helps to improve 
the x^-fit. 

For both baselines, the improvement gained by increasing the exposure is obvious to see - 
not only do the limits on sin^ 26'i3(true) and (5cp(true) decrease, as shown by the shrinking con- 
tours, but also the correlations between sin^ 26*13 (true) and (5cp(true) are decreased, as shown by 
the changing shapes of the contours from elliptical (for small sin^ 26*13 (true) or low exposure) to 
more circular. In the case that sin^ 2^13 (true) = 0.01, our first benchmark exposure will not be 
sufficient to place any reasonable bound on 5cp{ti\xe), although it can still produce a restriction 
on sin^ 26*13 (true). If sin^ 26*13 (true) = 0.1, the same exposure is able to define a 95% region of 
~ ±0.02 for sin^ 2^13 (true) and ~ ±45° for (5c'p(true). For our maximum exposure these regions 
shrink to ~ ±0.01 for sin^ 26^13 (true) and ~ ±20° for 5c'p(true). We will stress that these results 
will improve significantly if a prior on ^13 is included, as will all the results presented in this paper. 

7 Summary and conclusions 

Recent data from neutrino oscillation experiments indicate that the mixing angle ^13 is large, such 
that 0.01 < sin^ 2^13 < 0.1 at the 2a confidence level. Assuming that this hint will be confirmed 
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Total exposure 


(x lO^ipot-kt) 




Setup 


sin^ 26^13 (true) 


sin^ 2^13 (true) 


sin^ 2^^i3(true) 




0.01 


0.05 


0.1 




NH true 


IH true 


NH true 


IH true 


NH true 


IH true 


2290 km 


1590 (-) 


1811 (-) 


32 (80) 


50 (158) 


10 (25) 


20 (50) 


1540 km 


-(-) 


-(-) 


158 (502) 


250 (627) 


32 (100) 


50 (162) 



Table 5: Total exposure needed to achieve MH discovery with 100% coverage in (5cp(true) for CERN-Pyhasalmi 
and CERN-Slanic baselines. Results are presented for three different values of sin^ 2^13 (true) in the range of 0.01 
to 0.1. The numbers without (with) parentheses correspond to the 3ct (Scr) C.L. results. A dash indicates that 
this measurement is not possible with an exposure below our maximum limit of 3000 x lO^^pot-kt. 



sin^ 2^13 (true) 


Fraction of 5cp(true) at 3 a (%) 


Exposure 
300 X 10^1 pot-kt 


Exposure 
1000 X 10^1 pot-kt 


Exposure 
3000 X 10^1 pot-kt 


2290 km 


1540 km 


2290 km 


1540 km 


2290 km 


1540 km 


0.01 


0(0) 


0(0) 


26 (34) 


30 (31) 


59 (64) 


67 (66) 


0.05 


14 (31) 


36 (47) 


56 (63) 


65 (67) 


72 (74) 


75 (76) 


0.1 


26 (39) 


40 (46) 


59 (63) 


63 (64) 


71 (72) 


72 (73) 



Table 6: Comparison between CERN-Pyhasalmi and CERN-Slanic baselines for dcp fractions for which a 
discovery at 3ct C.L. is possible for CPV. Results are shown for three benchmark values of exposure and for three 
different choices of sin^ 20i3(true) in the range of 0.01 to 0.1. The numbers without (with) parantheses correspond 
to the results for a true NH (IH). 



with high statistical significance in the near future by further running of the T2K experiment, and 
by data from up-coming experiments, the optimization of next-generation long-baseline neutrino 
oscillation experiments must incorporate this knowledge. This is in contrast to previous studies 
which were performed without prior knowledge of ^13, and therefore had to take into consideration 
all possible values. 

In our study, we have focused on two complementary baselines: 2290 km, corresponding to 
the CERN to Pyhasalmi (Finland) baseline, and 1540 km, corresponding to the CERN to Slanic 
(Romania) baseline. This latter baseline is also close to the FNAL to DUSEL baseline (in the USA) 
of 1290 km. We have considered a CERN-based conventional wide-band beam which provides 



3 X 10^^ protons on target per year 30 , in combination with a Giant Liquid Argon detector. We 
did not fix the beam running time or the detector mass, in order to determine how the performances 
of both the setups vary with exposure (measured in units of pot-kt). We thus demonstrated that 
an incremental program, staged in terms of the exposure, can achieve the desired physics goals 
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within a realistically feasible timescale, whilst producing new, important and reliable results at 
each stage. Thereby the risks are minimized, as is the time between the start of the project and 
the point at which initial results appear. The program can be staged according to future technical 
developments and as a function of the discoveries achieved at each step, allowing the neutrino 
beam performance and the far detector mass to be developed in parallel. 

A relatively large value of ^13 and our choice of a long baseline ensure that it will be possible 
to determine the neutrino mass hierarchy with a 'pilot' Liquid Argon detector that has a mass 
in the range 10 to 20 kt, coupled to a conventional neutrino beam. The same setup will test at 
> 3 a" CL. if CPV is close to maximal and provide, in this case, a first measurement of the CP 
phase, 6cp- This approach is fully justified because at present, 6cp angles close to ±90° are not 
yet excluded, and these are the most easily accessible values. 

In the case that 6cp turns out to be inaccessible with the pilot setup, we will have gained the 
knowledge that Sqp lies in a region where we will need more ambitious exposures, achievable by 
either an increase of far detector mass and/or by a further increase of beam power. We have found 
that the reach corresponding to (50%) 70% coverage in 5c'p(true) for CPV will require an exposure 
of (1000)3000 X 10^^ pot-kt, ultimately limited by systematic errors. This 6cp{tTue) coverage is 
impressive when compared to ultimate scenarios, e.g. that of the multi-megawatt power Neutrino 
Factory which yields a coverage of 75% to 80% with a baseline magnetized detector of 100 kt. 

We found that the performance of the longer baseline, for mass hierarchy determination, is far 
superior to that of the shorter baseline. At 2290 km it is possible to obtain a 3 a measurement 
of the mass hierarchy, for all values of Scpitrue), for the full range 0.01 < sin^ 26*13 (true) ^ 0.1. 
The shorter baseline does not perform as well, especially for the smallest value of sin^ 26*13 (true) 
considered. We also find that for sin^ 26*13 (true) = 0.1, the shorter baseline requires roughly 3 
times more exposure than the longer baseline to make a 3 cr determination of the mass hierarchy. 
For a 5(7 measurement the difference is greater. If sin^ 2^13 (true) = 0.05 then the shorter baseline 
requires roughly 5 times as much exposure to make a 3 ex or 5 a measurement. Table |5] lists the 
precise exposures required by each of the baselines to achieve 100% coverage in (5cp (true) for mass 
hierarchy discovery, at the 3 a and 5 a confidence levels, for sin^ 26*13 (true) = 0.01, 0.05 and 0.1. 

In terms of CPV discovery potential, the shorter baseline requires relatively less exposure than 
the longer baseline to achieve the same 6cp{tTue) coverage, but the difference is not as great 
as for the hierarchy sensitivity and vanishes rapidly for the larger exposures and for growing 
sin^ 26*13 (true). To obtain a specified 5c'p(true)-fraction, the exposure needed for the 2290 km 
baseline is at most twice that required for the 1540 km baseline, but the difference is often much 
smaller. By defining experimental goals to be the discovery of CP violation at the 3 a confidence 
level, for 30%, 50% and 70% of the 6cp(tTue) parameter space for the current best-fit value of 
sin^ 26^13 (true) ~ 0.05, we defined three target exposures to define the stages of the experiment: 
an exposure of 300 x 10^^ pot-kt is roughly the exposure needed to obtain 30% coverage for both 
baselines and mass hierarchies, 1000 x 10^^ pot-kt gives 50% coverage and 3000 x 10^^ pot-kt 
provides roughly 70% coverage. Table |6] shows the CP coverage that can be obtained by each 
of the baselines for each of these benchmark exposure, at the 3 a and 5 a confidence levels, for 
sin^ 2^13 (true) = 0.01, 0.05 and 0.1. 
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Referring back to the hierarchy sensitivity in terms of these benchmark exposures, then with 
only our lowest exposure, the 2290 km baseline is able to make aba determination of the mass 
hierarchy for sin^ 26'i3(true) > 0.05 for all values of 6cp(tiue)] the 1540 km baseline can make a 
3cr measurement for these values of sin^ 2^13 (true). In other words, a 5 (3) cr result is guaranteed 
by the 2290 (1540) km baseline with our lowest exposure if sin^ 2^13 (true) > 0.05. Tables [s] and [6] 
summarize these results quantitatively. 

In this work we also quantitatively studied, for the first time, the optimization of the ratio of 
positive horn focusing (PHF) to negative horn focusing (NHF). Essentially, this is the same as the 
ratio of neutrino to anti-neutrino data. Firstly, we found that the optimal ratio for determining 
the mass hierarchy consists of roughly 30% PHF and 70% NHF, for both baselines, bearing in 
mind that we have to cover both the possibilities of a true normal hierarchy and a true inverted 
hierarchy. However, taking into account that NHF is experimentally disfavored, we suggest a 
mixture consisting of equal PHF and NHF for optimal sensitivity to the mass hierarchy. Once the 
mass hierarchy has been determined to a high confidence level, which will occur at an exposure 
well below that needed to obtain any significant CP sensitivity, the PHF to NHF ratio can then 
be altered, depending on the nature of the true hierarchy. A normal hierarchy favors more anti- 
neutrino data, so the 50 : 50 ratio should be maintained, whereas an inverted hierarchy prefers 
more neutrino data so that the NHF can be reduced to around 30% in this case. 

To conclude, we have found that if the current hints of large ^13 prove to be true, the measure- 
ment of the currently missing neutrino oscillation parameters - the neutrino mass ordering and 
the CP phase 6cp - will be achievable with an experimental program of smaller proportions and 
less ambitious design than those previously envisioned for the scenario of very small ^^13. To this 
end, we have demonstrated that a CERN-based conventional wide-band beam in conjunction with 
a Liquid Argon detector at the Pyhasalmi mine, at a baseline of 2290 km, would have exceptional 
sensitivity to the neutrino mass hierarchy and leptonic CP violation in the light of the emerging 
hints of large ^13. This setup fits well into an incremental program, staged in terms of the expo- 
sure, such that desired physics goals can be achieved within a realistically feasible timescale, whilst 
significant new results are produced at each stage. Staging the exposure is justified in view of 
the challenges involved in operating a high power superbeam and in building a giant underground 
neutrino detector - a staged approach is therefore effectively safer and possibly more cost-effective, 
whilst also being adaptable to future physics and technical developments. 
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